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Abstract 

A new paramctrization is introduced for the fixed point (FP) action 
in SU(3) gauge theory using fat links. We investigate its scaling 
properties by means of the static quark-antiquark potential and the 
dimensionless quantities roTc, Tcl\fo and r^^/a^ where Tc is the 
critical temperature of the deconfining phase transition, tq is the 
hadronic scale and a is the effective string tension. These quanti- 
ties scale even on lattices as coarse as a « 0.3 fm. We also mea- 
sure the glueball spectrum and obtain mo++ = 1627(83) MeV and 
m2++ = 2354(95) MeV for the masses of the scalar and tensor glue- 
balls, respectively. 
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1 Introduction 



One way to study quantum field theories beyond perturbation theory is to dis- 
cretize the Euclidean space-time, using the lattice spacing a as an ultraviolet 
regulator . Accordingly, the continuum action is replaced by some discretized 
lattice action. The basic assumption of universality means that the physical 
predictions - obtained in the continuum limit a —> - do not depend on the 
infinite variety of discretizing the action. At any finite lattice spacing, how- 
ever, the discretization introduces lattice artifacts. On dimensional grounds 
one expects that in purely bosonic theories these discretization errors go away 
as 0{a^), while in theories with fermions as 0{a). 

Sometimes it is stated that the lattice artifacts in pure Yang-Mills theories 
can be beaten simply by brute force - using the standard Wilson gauge action 
and a sufficiently small lattice spacing, i.e. by large computer power and memory. 
This is only partially true - for example when one calculates the pressure of a 
hot gluon plasma the computer cost grows like l/a^*^ therefore the size of lattice 
artifacts becomes crucial. 

Naturally, one can use the freedom in discretizing the action to minimize 
the artifacts. It has been shown by Symanzik || that the leading lattice 
artifacts can be cancelled in all orders of perturbation theory by tuning the 
coefficients of a few dimension d + 2 operators in bosonic theories (or dimension 
d -\- 1 operators for fermions). This improvement program can be extended to 
the non-perturbative regime 

A different approach, based on renormalization group (RG) ideas has 
been suggested in ref. 0. By solving the fixed point (FP) equations for asymp- 
totically free theories one obtains a classically perfect action ~ i.e. which has no 
lattice artifacts on the solutions of the lattice equations of motion. (One can 
say that the FP action is an on-shcU tree-level Symanzik improved action to all 
orders in a.) Although it is not quantum perfect, one expects the FP action 
to perform better in Monte Carlo simulations. This is indeed true in all cases 
investigated. 

The FP approach has been successfully applied to the two-dimensional non- 
finear cr-model @, || and the two-dimensional CP'^-model For SU(3) gauge 
theory the classically perfect FP action has been constructed and tested in 
[|o[ 0, [ij, |l^ and the ansatz has been extended to include FP actions for 
fermions as well 15, Jl6|. In the case of SU(2) gauge theory the FP action 
has been constructed in Jit] , [T^ , [l9t and its classical properties have been tested 
on classical instanton solutions, both in SU(2) and SU(3) pO[ . 

The FP action is not unique - it depends on the RG transformation chosen, 
and it is crucial to optimize the RG transformation to obtain an interaction 
range of the action as small as possible. The value of the FP action on a 
given field configuration can be calculated precisely by a classical saddle point 
equation. However, this step is too slow to embed into a Monte Carlo calculation 
and one has to invent a sufficiently fast but at the same time accurate enough 
method to calculate the FP action. 

In earlier works on SU(3) gauge theory the parametrization of the FP action 
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used Wilson loops and their powers. We investigate here a new, richer and more 
flexible parametrization using plaquettes of the original and of smeared ( "fat" ) 
links. This describes the FP action more accurately than the loop ansatz. For 
the RG transformation we choose the one investigated in ref. since the 
blocking kernel used there (and the resulting FP action) has better properties 
than the "standard" Swendsen blocking which uses long staples. Here we ap- 
proximate the same FP action with the new parametrization. We also improve 
the method of fixing the parameters: besides the known action values (for a 
given set of configurations) we also fit the known derivatives dA^^ {V)/SVfj,{n), 
i.e. we have a much larger set of constraints than previously. 

Although being much faster than the loop parametrization of comparable 
richness, this parametrization has a significant overhead compared to the Wilson 
action. Therefore it is not clear whether it is not better to use in pure gauge 
theory a faster but less accurate parametrization of the FP action. However, in 
QCD the cost associated with fermionic degrees of freedom dominates and one 
can afford a relatively expensive gauge action. In addition, with fermions it is 
much harder to decrease the lattice spacing, hence it could pay off to have a 
better gauge action as well. (Of course, because of the 0{a) artifacts, it is even 
more important to improve the fermionic part.) 

The paper is organized as follows. In section 2 we present the construction 
and parametrization of the FP action. Section 3 deals with the measurement 
of the critical couplings /3c of the deconfining phase transition corresponding to 
temporal extensions Nr = 2, 3 and 4, at various spatial volumes. In section 
4 we measure the static qq potential by using a correlation matrix between 
different (spatially) smeared gauge strings. In section 5 the scaling properties of 
the dimensionless quantities roTc, Tcj \fo and r^y/a are presented. In section 6 
the low lying glueball spectrum is measured in all symmetry channels. For the 
Wilson action the lowest lying 0++ state shows particularly large cut-off effects 
hence this quantity provides a non-trivial scaling test. Some technical details 
are collected in appendices A-D. 

2 A new parametrization of the FP action for 
SU(3) lattice gauge theory 

2.1 Introduction 

In this section we present a new ansatz for the parametrization which is very 
general and fiexible, and which allows to parametrize the FP action using more 
and more couplings without any further complications. The approach we use is 
building plaquettes from the original gauge links as well as from smeared ( "fat" ) 
links. In this manner we are able to reproduce the classical properties of the 
FP action better than with the loop parametrization. 

The new ansatz is motivated by the success of using fat links in simulations 
with fermionic Dirac operators ^ ^ |2j, |2^, ^ . Fat links are gauge links 
which are locally smeared over the lattice. In this way the unphysical short- 
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range fluctuations inherent in the gauge field configurations are averaged out 
and lattice artifacts are reduced dramatically. 

As mentioned above, earlier parametrizations of FP actions were based on 
powers of the traces of loop products along generic closed paths. Restricting 
the set of paths to loops up to length 8 which are fitting in a 2^ hypercube, 
one is still left with 28 topologically different loops some of them having 
a multiplicity as large as 384. In earlier production runs only loops up to 
length 6 (and their powers) have been used because including length 8 loops 
increases the computational cost by a factor of ^ 220 (l9|. Note that in the 
loop parametrization one needs length 8 loops to describe well small instanton 
solutions l^. 

The new parametrization presented here provides a way around these prob- 
lems, although the computational overhead is still considerable. We have calcu- 
lated the expense of the parametrized FP action and compared it to the expense 
of an optimized Wilson gauge code. The computational overhead amounts to a 
factor of ~ 60 per link update and comes mainly from recalculating the staples 
in the smeared links affected by the modified link. 



2.2 The FP action 

We consider SU(Af) pure gauge theory^ in four dimensional Euclidean space- 
time on a periodic lattice. The partition fimction is defined through 

Z(/3)= / dUe-^-^<-^\ (1) 



where dU is the invariant group measure and (3A{U) is some lattice regulariza- 
tion of the continuum action. We can perform a real space RG transformation, 



e 



^ / dU exp {-f3{A{U) + TiU, V))} , (2) 



where V is the blocked link variable and T(J7, V) is the blocking kernel defining 
the transformation, 

T{U,V) = E (ReTr(V;,(nB)Qt (ns)) - AA^^) . (3) 



N 

Here, Q^ijis) is a N x N matrix representing some mean of products of link 
variables U^{n) connecting the sites 2nB and 2{nB + jj) on the fine lattice 
and A/'^f is a normalization constant ensuring the invariance of the partition 
function. By optimizing the parameter k, it is possible to obtain an action on 
the coarse lattice which has a short interaction range. A simple choice for is 
the Swendsen blocking, which contains averaging over the 6 (long) staples along 
the direction /i. In rcf. JlSf the averaging was improved by including more paths 

^The following equations are given for general A^, although the numerical analysis and 
simulations are done for SU{3). 
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in Qp. The main idea of this block transformation is that, instead of using just 
simple staples, one additionally builds "diagonal staples" along the planar and 
spatial diagonal directions orthogonal to the link direction. In this way one 
achieves that each link on the fine lattice contributes to the averaging function 
and the block transformation represents a better averaging. In this paper we 
employ the RG transformation of ref. p^ , using, however, a completely new 
parametrization. 

On the critical surface at /3 — > oo, equation (Q) reduces to a saddle point 
problem representing an implicit equation for the FP action A^^ , 

^FP(y) = mm {^FP(C/) + nu, V)} . (4) 

The FP equation (kh can be studied analytically up to quadratic order in 
the vector potentials |[13|. However, for solving the FP equation on coarse 
configurations with large fluctuations one has to resort to numerical methods, 
and a sufficiently rich parametrization for the description of the solution is 
required. 



2.3 The parametrization 

In order to build a plaquette in the fiv-'plaiie from smeared links we introduce 
asymmetrically smeared links wif'' , where /i denotes the direction of the link 
and v specifies the plaquette-plane to which they are contributing. This asym- 
metric smearing suppresses staples which lie in the /ii^-plane relative to those in 
the orthogonal planes /xA, \^ v. Obviously, these two types of staples play a 
different role, and we know from the quadratic approximation and the numerical 
evaluation of the FP action that the interaction is concentrated strongly on the 
hypercube . Additional details on the smearing and the parametrization are 



given in appendix A.l 



From the asymmetrically smeared links we construct a "smeared plaquette 
variable" 

^z;^, =ReTr(l-W'Pi), (5) 
together with the ordinary Wilson plaquette variable 

= ReTr(l-C/Pi), (6) 

where 

W^^iH = W\^\n)W^>^\n + A)wMt(n + :>)Wi^)t(n), (7) 



and 

VfAn) = C/^(n)(7.(n + + v)Ul(n). (8) 
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Finally, the parametrized action has the form 
where we choose a polynomial in both plaquette variables, 

kl 

= PloU + poiW + p2ou'^ + PllUW + P02W'^ + . . . . (10) 

The coefficients pki, together with the parameters appearing in the smearing 



(cf. appendix A.l) should be chosen such that the resulting approximation to 
A^^ is sufficiently accurate. Note that the ansatz involves two types of param- 
eters: the coefficients in the asymmetric smearing enter non-linearly into the 
action while the coefficients pki enter linearly. 

For simulations with the FP action in physically interesting regions it is 
important to have a parametrization which is valid for gauge fields on coarse 
lattices, i.e. on typical rough configurations. We turn to this problem in the 
next section. 



2.4 The FP action on rough configurations 

The parametrization of the FP action on strongly fiuctuating fields is a difficult 
and delicate problem. In this section we describe briefly the procedure of obtain- 
ing a parametrization which uses only a compact set of parameters, but which 
describes the FP action still sufficiently well for the use in actual simulations. 
We also provide some details about the fitting procedure employed. 

In eq. (^) the minimizing field U = U{V) on the fine lattice is much smoother 
than the original field V on the coarse lattice - its action density is 30 — 40 
times smaller. This allows an iterative solution of eq. as follows. First one 
chooses a set of configurations V with sufficiently small fiuctuations such that 
for the corresponding fine field U{V) one can use on the rhs. an appropriate 
starting action (the Wilson action, or better the action ^o, eqs. ( A.15D and 



(A. 16) which describes well the FP action in the quadratic approximation). 
The values A^^{V) obtained this way have to be approximated sufficiently 
accurately by choosing the free parameters in the given ansatz. Once this is 
done, one considers a new, rougher set of configurations V for which the fields 
U(V) have fluctuations within the validity range of the previous parametrization 
and repeats the procedure described. After 3 such steps one reaches fluctuations 
typical for a lattice spacing o«0.2 — 0.3fm. 

In previous works only the action values A^^{V) have been fitted in order to 
optimize the parameters. Here we extend the set of requirements by including 
into the x^-function to be minimized the derivatives of the FP action with 
respect to the gauge links in a given colour direction a ~ 1, . . . , N'^ — 1, 

SA^^iV) 



(11) 
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Note that these derivatives are easily calculated from the minimizing configura- 
tion since they are given simply by ST{U, V)/SV^{n) at [/ = U{V). In this way 
one has 4(7V^ — 1)V local conditions for each configuration V instead of a single 
global condition, the total value A^^{V). (In addition, a good description of 
local changes is perhaps more relevant in a Monte Carlo simulation with local 
updates.) An important test for the fiexibility of the parametrization is whether 
both the requirements for fitting the derivatives and the action values can be 
met at the same time. This is indeed the case. 

For addressing questions concerning topology it is crucial that the parametri- 
zation describes accurately enough the exactly scale- invariant (lattice) instanton 
solutions jl^, |8[ |l^, 0. For this purpose we generate sets of SU(2) single 
instanton configurations embedded in SU(3) on a 12^ lattice with instanton 
radius p/a ranging from 3.0 down to 1.1. We then block the configurations 
down to a 6^ lattice (using the blocking which defines the RG transformation). 
As can be seen from figure 0, these are solutions of the FP equations of motion 
for radii p/a > 0.9.^ Note that for p/a < 0.9 the quantity T(t7, V) becomes non- 
zero, indicating that instantons of that size "fall through the lattice" , i.e. they 
are no longer solutions. The deviation from scaling at larger radii seen in figure 
|l| is due to the discontinuity at the boundary of the periodic lattice and is under 
control. 

In the final step, we first fit the derivatives on ^ 50 thermal configurations 
corresponding roughly to a Wilson critical coupling at « 2, « 5.1. In 
the following the non-linear parameters (defining the asymmetric smearing) are 
kept fixed, while we include in addition the action values and the derivatives of 
~ 75 thermal configurations at /3^^ = 2.8, 4.0, 7.0 and the action values of the 
instanton configurations. The corresponding is then minimized only in the 
linear parameters pki ■ 

To assure stability of the fit we check that the is stable on independent 
configurations which are not included in the fit. Using high order polynomials 
of the plaquette variables u and w there is a danger of generating fake valleys in 
the itw-plane (for u, w values which are not probed by the typical configurations 
and hence not restricted in x^)- The presence of such regions is dangerous since 
it can force the system to an atypical - e.g. antiferromagnetic - configuration. 
We find that this can be circumvented by choosing an appropriate set of pa- 
rameters. (Note that, as usually when parametrizing the FP action, there are 
flat directions in the parameter space along which the value of changes only 
slightly, i.e. there is a large freedom in choosing the actual parameters.) 

The smallest acceptable set of parameters consists of four non-linear pa- 
rameters describing the asymmetrically smeared links wji^'' and fourteen linear 
paramete rs pk i with < fc -I- / < 4. The values of these parameters are given in 



appendix A. 4 and fulfill the correct normalization. They form the final approx- 



imation of the FP action. 

*This is only approximately true: one should start from a very fine lattice and perform 
many blocking steps, and furthermore the periodic boundary conditions also violate (locally) 
the FP equations of motion. 
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Figure 1: Action values for SU(2) single instanton solutions on a G''^ lattice 
in units of the continuum instanton action Ainst as a function of the instanton 
radius p/a. Note, that A^^{V) = A^^{U) + T{U, V), where U is the minimized 
configuration on the fine lattice. 



Compared to the loop parametrization of ref. [[l3[ the present parametriza- 
tion gives a deviation from the true FP action values smaller by a factor of 2 
for configurations which are typical for the range of lattice spacing 0.03 fm < 
a < 0.2 fm. It also describes scale invariant lattice instantons for p/a « 1.1 to 
a precision better than 2%. Note that this parametrization is not intended to 
be used on extremely smooth configurations. In order to be able to describe 
the typical (large) fluctuations by a relatively simple ansatz we did not imple- 
ment the 0{a^) Symanzik conditions (cf. appendix |A.2| ) in the last step. In the 
intermediate steps (i.e. for smaller fluctuations), however, our parametrization 
(containing non-constant smearing coefficients ri{x), Ci{x)) is optimized under 
the constraint to satisfy the 0{a^) Symanzik conditions as well. 

In our Monte Carlo simulations we use only the final parametrization with 
constant rj and The parameters of the intermediate approximations to the 
FP action are not given here. 

To investigate the lattice artifacts with our parametrization, we perform a 
number of scaling tests, which are described in the subsequent sections. 
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3 The critical temperature of the deconfining 
phase transition 

3.1 Details of the simulation 

Using the parametrized FP action we perform a large number of simulations 
on lattices with temporal extension = 2, 3 and 4 at three to six different 
/3- values near the estimated critical /3c- Various spatial extensions N^/Nt = 
2.5. . .5 are explored with the intention of examining the finite size scaling of 
the critical couplings. Configurations are generated by alternating Metropolis 
and overrelaxation updates. 

In the equilibrated system we measure the Polyakov loops averaged over the 
whole lattice, 

^^i^ETr n (12) 

X t=0 

as well as the action value of the configuration after each sweep. Both values 
are stored for later use in a spectral density reweighting procedure. 

The details of the simulation and the run parameters are collected in tables 



C.l, C.2 and C.3, where we list the lattice size together with the values and 
the number of sweeps. However, near a phase transition the number of sweeps is 
an inadequate measure of the collected statistics, because the resulting error is 
strongly influenced by the persistence time and the critical slowing down. The 
persistence time of one phase is defined as the number of sweeps divided by the 
observed number of flip-flops between the two phases |^7|] . This quantity makes 
sense only for /3-values near the critical coupling (3c and has to be taken with 
care: for the small volumes which we explore, the fluctuations within one phase 
can be as large as the separation between the two phases, and the transition 
time from one state to the other is sometimes as large as the persistence time 
itself. The estimated persistence time Tp and the integrated autocorrelation 
time Tint of the Polyakov loop operator are listed in the last two columns in 



tables C.l, C.2 and C.2 



3.2 Details of the analysis 

For the determination of the critical couplings in the thermodynamic limit we 
resort to a two step procedure. First we determine the susceptibility of the 
Polyakov loops, 

XL^V,{{\m-{\L\r)^ K=iV3, (13) 

as a function of /3 for a given lattice size and locate the position of its maxi- 
mum. In the thermodynamic limit the susceptibility develops a delta function 
singularity at a first order phase transition. On a finite lattice the singularity is 
rounded off and the quantity reaches a peak value x^°^^ some PdVa)- 
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The critical coupling, i.e. the location of the susceptibility peak, is deter- 
mined by using the spectral density reweighting method, which enables the 
calculation of observables away from the values of /? at which the actual simu- 
lations are performed. This method has been first proposed in ^ and has 
been developed further by Ferrenberg and Swendsen ||3^, |3l| . 

In a second step we extrapolate the critical couplings for each value of Nr 
to infinite spatial volume using the finite size scaling law for a first order phase 
transition. 



(14) 



where /i « 0.1 is considered to be a universal quantity independent of Nj- p^ . 
In figure ^ we show the susceptibility xl as a function of f3. The solid lines 
are the interpolations obtained by the reweighting method and the dashed lines 
represent the bootstrap error band estimations. For the interpolations at a given 
lattice size we use the data at all beta values listed in tables C.l-C.S, although 
the runs at /3-values far away from the critical coupling do not infiuence the 
final result. 

In table |]we display the values of /3c(VJj) together with their extrapolations 
to VJj ^ oo according to formula ([l^). 





l3c{Nr - 2) 


MNr = 3) 


f3c{Nr = 4) 


6 


2.3552(24) 






8 


2.3585(12) 


2.6826(23) 




10 


2.3593(7) 


2.6816(12) 


2.9119(31) 


12 




2.6803(10) 


2.9173(20) 


14 






2.9222(20) 


00 


2.3606(13) 


2.6796(18) 


2.9273(35) 


h 


0.14(9) 


-0.05(7) 


0.25(9) 



Table 1: Results of the critical couplings /3c from the peak location of the 
Polyakov loop susceptibility and the corresponding infinite volume limit ob- 
tained according to relation (Q). The finite size scaling constant h is also 
given. 
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2.68 



2.69 



2.355 



2.360 
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Figure 2: The Polyakov loop susceptibility on lattices of size 4 x 12^, 3 x 12^, 
3 X 10^ and 2 x 10^. The solid curves are the interpolations using the spectral 
density reweighting, the dashed lines show the bootstrap error bands. 
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4 Scaling of the static quark-ant iquark potential 



4.1 Introduction 



An important part of lattice simulations is the determination of the actual 
lattice spacing a in order to convert dimensionless quantities measured on the 
lattice into physical units. From the static qq potential one usually determines 
a in units of the string tension a or in units of the hadronic scale tq. Using 
the string tension to set the scale is plagued by some difficulties. Since the 
noise/signal ratio increases rapidly with increasing r, the part of the potential 
V{r) from which the linear behaviour ar has to be extracted is measured with 
larger statistical errors. Further, due to the fact that the excited string has a 
small energy gap at large qq separations, it is difficult to resolve the ground 
state, and this can lead to a systematic error which increases the obtained value 
of the string tension. 

The use of the quantity rg circumvents these problems. In [ p3[ a hadronic 
scale Tc has been introduced through the force F{r) between static quarks at 
intermediate distances 0.2 fm ^ r < 1.0 fm, where one has best information 
available from phenomenological potential models jsj, Q and where one gets 
most reliable results on the lattice. One has 



(15) 



where originally |Q c — 1.65 has been chosen yielding a value rg « 0.49 fm = 
(395MeV)~^ from the potential models. However, on coarse lattices also this 
alternative way of setting the scale has its ambiguities as will be discussed in 



section 5.2 



Referring to precision measurements of the low-energy reference scale in 
quenched lattice QCD with the Wilson action 38 we collect values for 

c and Tc in table |^. The first line is calculated from data in |Q while the two 
last lines are taken from . 





c 


0.662(1) 


0.89 


1.00 


1.65 


1.65(1) 


4.00 


2.04(2) 


6.00 



Table 2: Parameter values for the determination of the hadronic scale through 
eq. (ITF 



The scaling of our parametrized FP action is examined by measuring the 
static qq potential and comparing the quantity tq {V{r) — y(ro)) versus r/rg at 
several values of (3. 

From the potential one can calculate rg and the effective string tension a. 
Finally one can test the scaling of the dimensionless combinations r^Tc, T^j \fa 
and r^y/a. This will be described in the next sections. 
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The scaling checks wih be pushed to the extreme by exploring the behaviour 
of the FP action on coarse configurations with very large fluctuations corre- 
sponding to a « 0.3 fm. This situation is not relevant for practical applications, 
and it becomes indeed more and more difficult to measure physical quantities 
due to the very small correlation length and rapidly vanishing signals. Never- 
theless, it is still interesting to investigate this situation in order to estimate 
the region in which the classical approximation to the renormalization group 
trajectory is still valid. 

For the standard plaquette action the static potential on coarse lattices shows 
strong violations of rotational symmetry and, before fitting it with 

a function of r, usually an empirical term (the lattice Coulomb potential minus 
1/r) is subtracted with an appropriate coefficient. On the contrary, for the 
FP action of ref. ||l^ , due to the proper choice of the RG transformation the 
resulting potential is (practically) rotational invariant]^. Here we do not aim at 
testing the rotational invariance of the potential but rather at determining ro 
and a. 



4.2 Details of the simulation 

We perform simulations with the FP action at six different /3-values, of which 
three correspond to the critical couplings determined in section ^. Configura- 
tions are updated by alternating Metropolis with overrelaxation sweeps. The 
spatial extent of the lattice s is c hosen to be at least ~ 1.5 fm, based on observa- 
tions in Table C^ contains the values of the couplings, the lattice 



sizes and the number of measurements. 

In order to enhance the overlap with the physical ground state of the po- 
tential we exploit smearing techniques. The smoothing of the spatial links has 
the effect of reducing excited-state contaminations in the correlation functions 
of the strings in the potential measurements. The operators which we measure 
in the simulations are constructed using the spatial smearing of . It consists 
of replacing every spatial link Uj{n), j = 1, 2, 3 by itself plus a sum of its 
neighbouring spatial staples and then projecting back to the nearest element in 
the SU(3) group: 



S^Uj{x) = Vsv(S)[Uj{x)+\sY.{Uk{x)Uj{x + k)ul{x+j) (16) 

+Ul{x ~k)Uj{x-k)Uk{x-k+j))y 



Here, 7^su(3)Q denotes the unique projection onto the SU(3) group element W 
which maximizes ReTr{WQ^) for an arbitrary 3x3 matrix Q. The smeared 
and SU(3) projected link SiUj{x) retains all the symmetry properties of the 
original link Uj{x) under gauge transformations, charge conjugation, refiections 

^Note that to cure the rotational invariance of the potential one has to improve not only 
the action but also the operators. For a more "rotational invariant" blocking the potential 
shows less violations of rotational symmetry. 
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and permutations of the coordinate axes. The set of spatially smeared hnks 
{SiUj{x)} forms the spatially smeared gauge field configuration. An operator 
O which is measured on a n-times iteratively smeared gauge field configuration 
is called an operator on smearing level Sn, or simply SnO. In the simulation 
of the static qq potential we use smearing levels with n = 0, 1, 2, 3, 4. The 
smearing parameter is chosen to be = 0.2 in all cases. 

The correlation matrix of spatially smeared strings is constructed in the 
following way. At fixed r we first form smeared string operators along the three 
spatial axes, connecting x with x + ri, 

SnVi{x,x + ri;T) = 

SnUi{x,T)SnUt{x + i,T) . . .SnUi{x + {r - l)i,T), 1=^1,2,3, (17) 

and unsmeared temporal links at fixed x, connecting r with t + t, 

Vi{T, T + t;x)^ U^ix, T)Ui{x, r + 1) . . . [^(f , t + {t ~ 1)). (18) 
Finally, the correlation matrix is given by 



Ci„i{r,t) = ( Ti-SiViix, X + ri; t)V4(t, t + t:x + ri) 

SraV^ix,x + rl; T + t)v}{T, T + t;x)Y (19) 

where (.) denotes the Monte Carlo average. In the following the correlation 
matrices are analyzed as described in section 4.3 



4.3 Details of the analysis and results 

In order to extract the physical scale through equation (|l^) we need an inter- 
polation of the potential and correspondingly the force between the quarks for 
arbitrary distances r. This interpolation of V{r) is achieved by fitting an ansatz 
of the form 

V{r) = + err (20) 

r 

to the measured potential values. 

We determine the scale in two steps. First we employ the variational tech- 
niques described in appendix ^ using the correlation matrix defined in eq. ( p^ 
for a given separation r. This method results in a linear combination of string 
operators SnV , n = 0, . . . ,4, which projects sufficiently well onto the ground 
state of the string, i.e. eliminates the closest excited string states. We then 
build a function using the covariance matrix which incorporates correlations 
between Cim{t) and Ci'm'{t')- Based on effective masses and on the val- 
ues we choose a region imin(?') < ^ < ^max('')- (Too small t values can distort 
the results due to higher states which arc not projected out sufficiently well. 
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while too large t values are useless due to large errors.) In this window we fit 
the ground state correlator by the exponential form Z{r) exp{—tV{r)), checking 
the stability under the variation of different parameters of the procedure. The 



results of these fits are collected in table C.6. 

A straightforward strategy is to fit the obtained values of V{r) by the ansatz 
(pO[). However, one can decrease the errors on a and a by exploiting the fact 
that the errors of V{r) at different values of r are correlated. Therefore, in the 
second step, we use the projectors to the ground state of the string for each r, 
obtained in the first step, and calculate the ground state correlator C{r, t) from 
the correlation matrix Cim.{r,t). After this we estimate the covariance matrix 
Cov(r, t;r',t') from the bootstrap samples of C{r,t) and use this to build a 
function, fitting C{r,t) with the expression Z(r) exp(— t(Vb — a/r + err)). We 
use the fit range tminir) < t < tmaxir) determined in the first step. The fit 
range in r is determined by examining the values and the stability of the 



fitted parameters. The results of these fits are given in table |C.5 , 

Having in hand a global interpolation of the static potential for each /3-value, 
we are able to determine the hadronic scale rg in units of the lattice spacing 
through eq. (|l^). The value of c is chosen appropriate to the coarseness of the 
lattice and the fit range in r. 

In addition, we repeat the second step, but restricting this time the values 
of r to values close to to have a local fit to V{r) |8j. This fit is used then 
again to determine Vc/a (and accordingly tq/o) from the relation (p^. 

The final results for r^/a are listed in table ^ where the first error denotes the 
purely statistical error. The second one represents an estimate of the systematic 
error and marks the minimal and maximal value of tq /a obtained with different 
local fit ranges and different reasonably chosen values of c. These ambiguities 



are discussed in detail in section 5.2 



p 




ro/a 


3.400 




4.833(39)(+22) 


3.150 




3.717(23)(+^i^) 


2.927 


4 


2.969(14) (+f J 


2.860 




2.740(10)(+3[) 


2.680 


3 


2.237(7)(+J^) 


2.361 


2 


1.500(5)(+?^) 



Table 3: The hadronic scale r^/a determined from local fits to the potential. 
The first error denotes the statistical error and the second is the estimate of the 
systematic error. 



In figure y we display the potential values. The dashed line is obtained by a 
simultaneous fit to all the data respecting the previously chosen fit ranges in r. 
The dotted line representing the result of |4^ (obtained on anisotropic lattices 



14 



using an improved action) is hardly distinguishable from our dashed line. 
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Figure 3: Scaling of the static qq potential V(r) expressed in terms of the 
hadronic scale tq. The unphysical constant roViro) is subtracted for each lattice 
spacing such that the curves match at r/ro = 1. The dashed line is a fit of the 
form ( pO|) t o the data. The dotted line, which practically coincides with our fit, 
is from ||43[. 

It could be useful to have an empirical interpolating formula connecting the 
lattice spacing to the bare coupling. (Analogous fits for the Wilson action are 
given in refs. @-) The expression 

ln(a/ro) = -1.1622(24) - 1.0848(95)(/3 - 3) + 0.156(17)(/3 - 3)^ (21) 

describes well the data points in the range 2.361 < (3 < 3.4. The fit is shown in 
figure H 

5 Scaling of the critical temperature and ro\/^ 

To further study the scaling properties of the FP action we examine the dimen- 
sionless combinations of physical quantities Td \fd-, r^Tc and r^^fa. 

In this section we present and discuss the results for the FP action and 
compare them to results obtained for the Wilson action and different improved 
actions whenever it is possible. 
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Figure 4: The measured data points of ln(a/ro) (circles) and their phenomeno- 
logical description in terms of a polynomial quadratic in (3 (solid line). The 
plotted points are the values of ro/a from local fits. 



5.1 Te/V^ 

Let us first look at the ratio 
the string tension]^. 



fj, the deconfining temperature in terms of 



action 



FP action 



Wilson ^2 1 

1 X 2 

DBW2 " 
Iwasaki 
Bliss 



2.927 0.624(7) 
2.680 0.622(8) 
2.361 0.628(11) 




oo 
oo 
oo 
oo 
oo 



0.630(5) 

0.634(8) 

0.627(12) 

0.651(12) 

0.659(8) 



Table 4: Results of the deconfining temperature in units of the string tension 
obtained with the FP action and continuum values for different other actions. 



In table ^ we collect all available continuum extrapolations together with 
the results for the FP action. The data obtained with the Wilson action is 
taken from where they use the Tc values at Nr = 4 and 6 from and 
extrapolate finite volume data for Tc at Nr — 8 and 12 from |^ to infinite 
volume. For the value of a they use the string tension parametrization given 
in [p^. The data for the 1x2 tree level improved action is again taken from 



^Here and in the following we refer to the quantity cr obtained from the three parameter 
fits of the form (|20|) to the static potential as the (effective) string tension. 
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[ p2| . The data denoted by RG improved action is obtained with the Iwasaki 
action ^] and is taken from Q. We also include the value by Bliss et al. [Q 
from a tree level and tadpole improved action. Finally we quote the results 
from the QCD-TARO collaboration @ obtained with the DBW2 action|]. The 
extrapolations to the continuum stem from p9| where a careful reanalysis has 
been done. 

For extracting the string tension we follow a simple approach. As described 
above, we perform fits to the on-axis potential values only and therefore we are 
limited to a small number of different fitting ranges. Nevertheless, the values of 
a obtained this way and quoted in table C.5 are stable and vary only within their 
statistical errors over the sets of sensibly considered fit ranges. However, the 
error on tr changes considerably, i.e. up to a factor of 5, depending on whether 
distance r = 1 is taken into account or not. Just to play safe we neglect distance 
r = 1 in the fits, even if the would allow it. 

The values are displayed in figure ^ together with the data as mentioned 
above. Our data is compatible within one standard deviation with the contin- 
uum extrapolation of the Wilson data and we observe scaling of the FP action 
within the statistical errors over the whole range of coarse lattices corresponding 
to values of = 2, 3 and 4. 



0.70 



0.65 



<3<]<i B 

rn -J __ 



0.60 - 



O RG improved action 

□ DBW2 action 

1x2 tree level improved 

■O Wilson action 

■ FP action 



0.00 0.05 0.10 0.15 



0.20 0.25 



0.30 



Figure 5: T^j \fo vs. for different actions. 



^DBW2 means "doubly blocked from Wilson in two coupling space" 
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5.2 roTe 



Unfortunately, precise determinations of rg/a are missing in the literature except 
for the Wilson action Q and, in contrast to Tc/\/d, we are not able to 
compare our data to other actions such as the Iwasaki, DBW2 or the 1x2 
tree level improved action. In fact, the determination of rg/a is a delicate issue 
and systematic effects due to different methods of calculating the force can be 
sizeable. Due to the fact that extracting the derivative of the potential from 
a discrete set of points is not unique, the intrinsic systematic uncertainty is 
not negligible at intermediate and coarse lattice spacings a > 0.15 fm. For 
example, in an accurate scale determination of the Wilson gauge action in pTf 
the authors quote a value of ro/a — 2.990(24) at I3w = 5.7. This is to be 
compared with r^/a = 2.922(9) of ref. for the same action and /3-value. In 
view of the claim in to have included all systematic errors and the high 
relative accuracy (~ 0.3%) of the data in |Q, this systematic difference seems 
to be a serious discrepancy. Even on finer lattices there are ambiguities: at 
(3w — 6.2 the authors of |Q obtain r^/a = 7.38(3), while in a value of 
ro/a = 7.29(4) is quoted. 

In that sense our results concerning rpTc have to be taken with appropriate 
care. In table || we collect the data for tqTc from our measurements with the FP 
action together with the data from measurements with the Wilson action. The 
critical couplings corresponding to Nj- = 4, 6, 8 and 12 are taken from js^ while 
the values for ro/a are from the interpolating formula in [3^ . The quoted errors 
are purely statistical. The continuum value is our own extrapolation obtained 
by performing a fit linear in the leading correction term 1 /iV^ and discarding the 
data point at Nt — 4. Finally, the values are plotted in figure |6| for comparison. 



Nr 


Wilson action 


FP action 


2 




0.750(3) 


3 




0.746(3) 


4 


0.719(2) 


0.742(4) 


6 


0.739(3) 




8 


0.745(3) 




12 


0.746(4) 




oo 


0.750(5) 





Table 5: Results for the critical temperature in terms of the hadronic scale, 
roTc, from measurements with the Wilson action and the FP action. 



The Wilson action shows scaling violation for roTc of about 4% at = 4, 
while at = 6 it is already smaller than about 1.5%. In that sense this quantity 
provides a high precision scaling test and thus a very accurate computation of 
the low-energy reference scale ro/a on the 0.5% level is of crucial importance. 
The lack of data for different actions is an indication that this is indeed a difficult 
task. Although the required statistics is in principle accessible to us, we do not 



18 



0.77 



0.76 - 



0.75 - Ot 



O 



o 



0.74 



0.73 - 



0.72 - 



0.71 



O 



0.00 0.05 



0.10 0.15 
1/N' 



O Wilson action 
■ FP action 



0.20 0.25 0.30 



Figure 6: tqTc vs. for the Wilson and the FP action. The empty circles 

represent data from measurements with the Wilson action and the filled squares 
denote the results obtained with the FP action. The solid error bars show the 
purely statistical error, while the dashed ones indicate the systematic error from 
the ambiguities in determining the force on coarse lattices. 
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have full control over the systematic ambiguities in the calculation oir^/a on the 
required accuracy level. Nevertheless we observe in principle excellent scaling 
within 1% or two standard deviations for the FP action even on coarse lattices 
corresponding to Nt = 3 and 2, however, this statement is moderated in view 
of the large systematic uncertainties. 



5.3 ro\/CT 

To obtain the dimensionlcss product ro \/a we use the values of ro /a in table |^ 
obtained from local fits and the values of a as determined in section 5.1, where 
cr is determined from the long range properties of the potential. 

In table || we collect the resulting values of tq^/u. We can extrapolate 
to the continuum by performing a fit linear in {a/ro)^ and obtain r^^/a = 
1.193(10). For comparison we calculate the data for the Wilson action from 
the interpolating formula for ro/a in and the string tension parametrization 
in ||3^. The continuum extrapolation for the Wilson data is taken from the 
analysis of Teper in M] . 



Wilson action 


FP action 


5.6925 1.148(12) 
5.8941 1.170(19) 
6.0624 1.183(13) 
6.3380 1.185(11) 

oo 1.197(11) 


2.361 1.194(21) 
2.680 1.196(15) 
2.860 1.190(23) 
2.927 1.191(12) 
3.150 1.185(16) 
3.400 1.198(12) 
oo 1.193(10) 



Table 6: rov^ for the Wilson and the FP action. 



Figure shows the scaling behaviour of for the Wilson action (empty 

circles) and the FP action (filled squares) as a function of (a/ro)^. The error 
bars are purely statistical and are dominated by the uncertainty from the string 
tension. Therefore the systematic ambiguities present in ro/a are not visible 
within the error bars. 

The Wilson action shows a scaling violation of about 4% ai (3 — 5.6925 
{Nt = 4), while no scaling violation is seen for the FP action even on lattices 
as coarse as (3 — 2.361 {N^ = 2). We would like to emphasize that this is a 
non-trivial result, since ro/a and ^/aa are determined independently of each 
other. However, with the data presently available to us it is difficult to extract 
the string tension with the accuracy needed to see a striking difference to the 
Wilson action for /3- values corresponding to Nr > 4. This is mainly due to the 
lack of measurements of the off-diagonal potential values. 
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Figure 7: Scaling behaviour of r^yju for the Wilson action (empty circles) and 
the FP action (filled squares). 
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6 Glueballs 



6.1 Introduction 

Glueballs are the one-particle states of SU(iV) gauge theory. They are charac- 
terized by the quantum numbers J^'^, denoting the symmetry properties with 
respect to the 0(3) rotations (spin), spatial reflection and charge conjugation. 
However, the lattice regularization does not preserve the continuous 0(3) sym- 
metry, only its discrete cubic subgroup, therefore the eigenstates of the transfer 
matrix are classified according to irreducible representations of the cubic group. 
There are five such representations: Ai, A2, E, Ti, T2, of dimensions 1, 1, 2, 
3, 3, respectively. Their transformation properties can be described by poly- 
nomials in X, y, z as follows: Ai ~ {1}, A2 ~ {xyz}, E ~ {x^ — — z^}, 
T\ ^ {a;,j/,z} and ~ {xy^xz^yz^^ where a;, y, z are components of an 0(3) 
vector. In general, an 0(3) representation with spin J splits into several repre- 
sentations of the cubic group. Looking at the corresponding polynomials, it is 
rather obvious that the splitting starts at J = 2: (J = 0) ^ Ai, (J = 1) ^ Ti, 
( J = 2) ^ {E,T2). The full 0(3) rotation symmetry is expected to be restored 
in the continuum limit. This restoration manifests itself e.g. in the fact that a 
doublet E and a triplet T2 (for a given choice of quantum numbers PC) become 
degenerate to form together the J = 2 states with 5 possible polarizations. 

The main obstacle in the computations of glueball masses on the lattice is 
the fast decay of the signal in the correlation functions of the gluonic excitations, 
due to the fact that the glueball masses are relatively large (mg ^1-6 GeV). 
For this reason a small lattice spacing a is required to follow the signal long 
enough. On the other hand, the physical lattice volume should be larger than 
L > 1.2 fm to avoid finite size effects. This finally results in a large L/a 
making it hard to obtain the statistics which is usually required. One possible 
way around this dilemma is the use of anisotropic lattice actions, which have 
a finer resolution in time direction, <ti aa-, and where one can follow the 
signal over a larger number of time slices. Although this idea is not new | pT[ , 
it has been revived only recently by Morningstar and Peardon ||5^, Using 
an anisotropic improved lattice action they investigated the glueball spectrum 
below 4 GeV in pure SU(3) gauge theory and improved the determinations of the 
glueball masses considerably compared to previous Wilson action calculations. 
Recent calculations with the Wilson action comprehend works by the UKQCD 
collaboration Q and the GFll group It can be said that all three 

calculations are in reasonable agreement on the masses of the two lowest lying 
0++ and 2++ glueballs. 

Despite this agreement, Wilson action calculations of the 0++ glueball mass 
show huge lattice artifacts of around 40 % at a w 0.15 fm and still 20 % even 
at a « 0.10 fm. From this point of view the O"*""*" glueball mass is particularly 
interesting, besides its physical relevance, since it provides an excellent test 
object on which the scaling behaviour of different actions can be checked and 
the achieved reduction of discretization errors can be sized. In this sense let us 
emphasize that our intention here is twofold: firstly, our calculation provides a 
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new and independent determination of glueball masses using FP actions, and 
secondly, we aim at using the glueball spectrum, in particular the mass of the 
0++ glueball, as another scaling test of the FP action. Although we observe 
that the FP action scales well in quantities like rgTc, Td \fa or r^y/a^ lattice 
artifacts could be, in principle, quite different for other physical quantities, in 
particular rouiG or mGl\fo. 



This section is organized as follows. In subsection 6.2 we describe the details 
of the simulations including the generation of the gauge field configurations and 
the measurements of the operators. The extraction of masses from the Monte 



Carlo estimates of glue ball correlation functions is described in subsection p. 3 . 
Finally, subsection |6.4| contains the results of our glueball measurements. 

6.2 Details of the simulation 

We perform simulations at three different lattice spacings in the range 0.1 fm < 
a < 0.18 fm and volumes between (1.4 f m)'^ and (1.8 fm)'^. The simulation 



parameters for our runs are given in table C.7 . 

The gauge field configurations are updated by performing compound sweeps 
consisting of alternating over-relaxation and standard Metropolis sweeps. 

First, a rather small preliminary simulation at (3 ~ 2.86 is performed. Using 
the results of some pilot runs, we determine a set of five loop shapes which have 
large contributions to the A^^ channel. Using the labelling of Berg and Billoire 
\ ^t\ these are the length-8 loop shapes 2, 4, 7, 10, 18. They are measured on 
five smearing levels 5„, n = 2, 4, . . . ,10 with smearing parametei]^ — 0.2 and 
subsequently projected into the Af^ channel. 

In the two large simulations a,t (3 — 3.15 and 3.40 we measure all 22 Wilson 
loop shapes up to length eight (see ) on the same smearing levels mentioned 
before and project them into all 20 irreducible glueball channels. 

A considerable part of the simulation time is used to measure all the 22 loop 
shapes. Some of them may turn out to be superfluous in the sense that they 
give a much worse signal/noise ratio than the others. On the other hand, one is 
interested in having a set of operators as large as possible to build up the wave 
function of the lowest glueball state (more precisely, to cancel the unwanted 
contributions from the neighbouring states in the spectrum) . Having measured 
all these operators will allow us to identify the important loop shapes to be used 
in future simulations. 

The projections of the loop shapes into the irreducible representations of 
the cubic group are done according to the descriptions in [ p7| , [5^ . The corre- 
lation matrix elements are then constructed from the projected operators and 
Monte Carlo estimates are obtained by averaging the measurements in each 
bin. We measure all possible polarization^ in a given channel and add them 
together in the correlation matrix. This eventually suppresses the statistical 

*For details of the smearing we refer to subsection |4.2| . 

^In analogy to choosing different magnetic quantum numbers m for given angular momen- 
tum I in the 0(3) group. 
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noise more than just increasing the statistics since the different polarizations 
are anti-correlated. 

For the extraction of the glueball masses in the A^^ representation (which 
has the same quantum numbers as the vacuum) one has to consider vacuum- 
subtracted operators. For this purpose we also measure the expectation values 
of all the operators. 



6.3 Details of the analysis 

The glueball masses are extracted using the variational techniques described in 
appendix Let us put some remarks which are related to the analysis of the 
glueball masses in particular. 

As we are measuring a large number of operators (up to 145), normally 
some of them contain large statistical noise. Therefore we only keep a set of 
well measured operators, on which the whole procedure is numerically stable 
and well defined. 

Another remark concerns the vacuum subtraction necessary in the Af^ 
channel. To obtain vacuum-subtracted operators one usually considers 
(jf^^{T) — (f){T) — (O|0(r)|O). However, we follow a different strategy and treat 
the vacuum on the same footing as the other states in the vacuum channel. As 
it turns out, the vacuum state can be separated in this procedure with very 
high accuracy and it is safe to consider only the operator basis orthogonal to 
the vacuum in the fitting procedure. For this purpose we cut out the vacuum 



state obtained from solving the generalized eigenvalue equation (B.2), i.e. we 
only consider the correlation matrix^ 

Cf,{t)^{v,,C^'{t)v,), (22) 

with I, j running from i, j ~ 2, . . . ^ K < M in the further analysis (« = 1 being 
the vacuum state). In our experience this strategy yields the most stable sub- 
traction of the vacuum contribution with respect to the statistical fluctuations 
of the subtracted operators. 

In the last step for extracting the glueball masses the large correlation matrix 
is truncated toalxlor2x2 matrix, which is subsequently fitted in the fit 
range ^min • • ■ imax taking both temporal correlations and correlations among 
the operators into account. The corresponding covariance matrix is calculated 
from jackknife samples and the error is estimated using a jackknife procedure. 
The choice of tmax is not crucial and is usually taken according to the relative 
error of the matrix elements under consideration and the x^-function. More 
important is the correct choice of tmin- Since excited glueball states are rather 
heavy we do not expect large contamination of the ground state correlators from 
excited states even on time slice t = 1 and therefore tmin = 1 is usually chosen. 
In particular this choice is safe if we fix = 1 and ti = 2 rather than to = 



and ti = 1 in the variational method, eq. (B.3). Indeed, in the former case the 
X^-function remains more stable when we increase ^min 

^"See appendix |B|for notations. 
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check for the consistency of the resulting masses (as an example take the results 
in table Icj). 



6.4 Results 



The results of the fits to the glueball correlators are collected in the appendix 
in tables C.S - C.IC. We include the results of different fitting ranges in t in 
the tables in order to give an impression of the stability of the fits. In each 
channel the result highlighted in boldface is our final choice and represents a 
most reasonable mass for the given channel. These final mass estimates in units 



of the lattice spacing are collected in table C.ll. 

To compare the values it is convenient to use tq to set the scale. In table C.12| 



we list our estimates of the glueball masses expressed in terms of tq , while figure 
^ and ^ show our values for the A^^ and the E'^'^ , T^'^ channels, respectively, 
together with results from different calculations with the Wilson action (crosses) 
[I9I p3, p6| and the calculations of Morningstar and Peardon ||5^, 53 and Liu 



[ jQI with a tree level/tadpole improved anisotropic action (empty symbols). 
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Figure 8: Glueball mass estimates for the A^'^ channel. Results from simula- 
tions of the Wilson action (crosses) and a tree level/tadpole improved anisotropic 
action (empty circles) are shown together with the results obtained with the FP 
action (filled circles). 
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To compare our results to the continuum values of the various collaborations 
we resort to where the above Wilson action results have been expressed in or 
converted to units of rg using the interpolating formula for the Wilson action |Q 
and, whenever necessary, the continuum extrapolation has been redone. Our 
continuum result for the 0++ glueball mass is an extrapolation to the continuum 
using a fit function linear in (a/ro)^. The data in the other channels does not 
allow to do an extrapolation, thus we simply quote the masses obtained on the 
finest lattice (a = 0.10 fm) in brackets. The comparison of our results to the 
continuum values of the other groups is listed in tables and ||. 

Note that one observes restoration of the degeneracy within the statistical 
errors for the 2++ state as well as for the 2 state. All our mass estimates 
agree with the best earlier results within the statistical errors. 



CoUab. 




romo++ 


rom2++ 


year 


UKQC 
Teper 
GFll 1 
M&P 1 
Liu [59 


D 
49 
56 
53 


54 4.05(16) 
4.35(11) 
4.33(10) 
4.21(11)(4) 
4.23(22) 


5.84(18) 

6.18(21) 

6.04(18) 

5.85(2)(6) 

5.85(23) 


1993 
1998 
1999 
1999 
2000 


FP action 4.12(21) 


[5.96(24)] 


2000 



Table 7: Comparison of the two lowest glueball masses in units of tq. Our 2+"'" 
value is not extrapolated to the continuum but is the mass obtained at a lattice 
spacing a = 0.10 fm. 



CoUab. 




romo-+ 


rom2-+ 


romi+- 


year 


Teper p9 




5.94(68) 


8.42(78) 


7.84(62) 


1998 


M&P [53 




6.33(7)(6) 


7.55(3) (8) 


7.18(4)(7) 


1999 


FP action 


[6.74(42)] 


[8.00(35)] 


[7.93(78)] 


2000 



Table 8: Comparison of glueball masses in units of tq. Values in brackets denote 
masses obtained at a lattice spacing a = 0.10 fm and are not extrapolated to 
the continuum. 



Finally we convert the scalar and tensor glueball masses into physical units 
using ro ~ 0.49 fm = (395MeV)-i. We obtain 1627(83) MeV for the 0++ and 
2354(95) MeV for the 2++ glueball mass, respectively. Note that the latter value 
corresponds to the glueball mass measured at a lattice spacing a = 0.10 fm. 

As mentioned in the introduction, it is well known that glueball masses are 
difficult to measure on the lattice. Indeed, we can barely resolve higher lying 
glueball states and measuring excited states becomes impossible at the lattice 
spacings currently available to us. In this sense we can not really take advantage 
of the parametrized FP action, which is intended to be used on coarse lattices. 
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Figure 9: Glueball mass estimates for the 2++ channel. Results from simulations 
of the Wilson action (crosses) and a tree level/tadpole improved anisotropic 

action (empty symbols) arc shown together with the results obtained with the 
FP action (filled symbols). Squares and circles denote the E'^'^ and mass 
estimates, respectively. 
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One way around this difficulty is the use of anisotropic lattices, where the 
lattice spacing in temporal direction is much smaller than in spatial direction, 
Or <C Co- • The work on the application of the FP approach to anisotropic lattice 
gauge actions is in progress pl[ |. 

7 Summary 

In this work we have presented a new parametrization of the FP action of a 
specific RGT. It uses simple plaquettes built from single gauge links as well as 
from smeared ("fat") links. It reproduces the classical properties of the action 
excellently and respects approximate scale invariance of instanton solutions. 
Since in addition to the FP action values we parametrize the derivatives with 
respect to the gauge fields, local changes of the action in a MC simulation are 
better represented. 

The parametrization has been optimized at lattice spacings suitable for per- 
forming simulations on coarse lattices up to a « 0.3 fm. 

For subjecting the action to scaling tests we have determined its critical 
couplings /3c on lattices with temporal extensions Nr = 2, 3 and 4. For each 
Nr we have performed simulations on several lattices for a finite size scaling 
study. Furthermore, we have measured the static quark-antiquark potential 
at various values of the gauge coupling corresponding to a ~ 0.1 — 0.3 fm. 
From the potential we have extracted the commonly used reference scale tq 
and the effective string tension a in order to check the scaling behaviour of the 
parametrized FP action by means of the dimensionless quantities ^oTcj T^j ^fa 
and TQ^fa. In all the quantities we observe excellent scaling within the statistical 
errors, even on our coarsest lattices. 

Additionally, we have measured the glueball spectrum in all symmetry chan- 
nels. The A^'^ channel, which shows particularly large lattice artifacts in mea- 
surements with the Wilson gauge action, is an excellent candidate for testing the 
improvements achieved with the parametrized FP action. We observe scaling 
of the glueball masses and restoration of the rotational symmetry in the 2++ 
and 2 ^ channel within the statistical errors. For the glueball masses we obtain 
1627(83) MeV for the 0++ glueball in the continuum and 2354(95) MeV for the 
2++ glueball at a lattice spacing of a = 0.10 fm. 

Acknovifledgements We would like to thank Peter Hasenfratz for useful sug- 
gestions and discussions. 
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A The parametrization 

A.l Details of the parametrization 

Let us introduce the notation sj^''' (n) for the sum of two staples of gauge hnks 
in direction /i in the /xi^-plane: 

+ Ul{n - i))U^{n - v)U,{n - + fi) . (A.l) 

Besides the usual symmetric smearing, we shall also use a non-symmetric 
smearing. For the symmetric smearing define 



rd0 



anc^ 

x^{n)=ReTr{Ql{n)Ul{n)). (A.3) 

To build a plaquette in the /nj^-planc from smeared links wc introduce asym- 
metrically smeared links. First definepl 



^,)]U,. (A.4) 



Using these matrices we build the asymmetrically smeared links 

=U, + ci(x^)QM + c2(x^)QMt/tQM + .... (A.5) 
Here 77(x), Ci{x) are polynomials: 

?7(x^) = r/(") + 77(1)3;^ + r,(2)x2 + . . . (A.6) 

and 

c,ix^ = c^^+^^x,+c^^xl + .... (A.7) 
A.2 The C(a^) Symanzik conditions 

In this appendix we derive the O(a^) Symanzik conditions [|[ ^ ^ ^ ^ 
by considering constant non-Abelian gauge potentials. The formulas apply to 
general SU(A^). 

^^Note that a;^(n) is negative: —4.5 < x^(n) < 0. 
^^The argument n is suppressed in the following. 
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In the continuum we have one scalar, gauge-invariant dimension-4 operator 
^o = -^ETr(.^^.), (A.8) 

and three dimension-6 operators: 

Ri = \Y.^i-{{D,J^^,f), (A.9) 

R2 = lY.^i-{{D^J'.xf), (A.IO) 

Rs^IY. {Df^T^xD^T^x) ■ (A.ll) 

According to Symanzik H, ||] the 0{a^) lattice artifacts of an action are 
described by an effective continuum action where to the usual continuum action 
(oc Rq) additional terms proportional to a^Ri, a^R2 and a^Rs are added with 
appropriate coefficients. In fact, the equations of motion are D^T^x = 
hence the term with R^ can be eliminated by a change of variables, hence does 
not affect the lattice artifacts in on-shell quantities, e.g. masses. (Note that 
the static qq potential is an off-shell quantity, it depends on the choice of the 
operators. For such quantities one has to improve the operators as well to get rid 
of artifacts.) The 0{a^) lattice artifacts of on-shell quantities can be eliminated 
in all orders of perturbation theory by adding to the original lattice action two 
additional terms which in the naive continuum limit are proportional to a^Ri 
and a^i?2, with appropriate coefficients. On the tree level the absence of 0{a?') 
artifacts means that when one expands the lattice action in powers of a, for 
smooth fields the coefficients of i?i and R2 vanish. The coefficient of R^ is not 
required to vanish (and usually it does not for the FP action). 



For the specific lattice gauge action ansatz considered in section 2.3 one 
obtain^ 



1 



= |i?o(l + (4 + 277)ci) 

IJ,<v 



+ ^Ri (1 - 2ci(l - iv) + ^(1 - - 2c2) 

+ \r3 (ci + i(l + 2ij){cl - 2c2)) . (A.12) 

The normalization condition is obtained from the coefficient of Rq^ 

Pio+Poi(l + (4 + 2r7)ci) = 1. (A.13) 



^''Prom the non-linear parameters only the zeroth order coefficients contribute to the normal- 
ition and the 0{a^) Symanzik condition, 
and jyf") — > 77 in the rest of this subsection. 



ization and the 0(a^) Symanzik condition. To keep notation simple we substitute c'"-' 
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The first C(a^) Symanzik condition requires the coefficient of Ri to vanish, 

Pio+Poi (^l~2ci{l-A7j) + ^{l-rjf{cl-2c2)^ =0. (A.14) 

It is interesting to see that the operator R2 is absent and hence the second 0(0^) 
Symanzik condition is satisfied automatically for the general ansatz considered 
here. (Note that when the FP action is expressed in terms of simple loops some 
of them give a nonzero contribution to R2I) 



A. 3 The quadratic approximation 

The couplings of the FP action can be calculated analytically in the quadratic 
approximation ||ll|, |l^. By fitting the leading order nonlinear parameters r/'"', 

cf^ , £^2^^ and pio , poi to the quadratic part of the FP action we can check the 
flexibility and the quality of the parametrization. Although the true FP action 
fulfllls the tree-level Symanzik conditions to all orders in a, an approximate 
parametrization introduces small violations of all these conditions. However, 
one can exploit the freedom in the parametrization to correct for this and to 
fulflU explicitly the 0{a^) on-shcll Symanzik conditions. The linear parameters 
Pio and poi are determined as functions of 77'*'-', cf^ and Cj^'' by eqs. (A. 12) and 



(A.14). The fit to the exactly known quadratic approximation of the FP action 



yields the following result for the three nonlinear parameters: 

ryf") = 0.082 , c^°^ = 0.282 , = 0.054 , (A.15) 

with the corresponding plaquette coefficients 

Pio = -0.368095 , poi = 0.629227 . (A.16) 

This action is denoted by Aq and is a good approximation to the FP action for 
sufficiently smooth fields. 



A. 4 The parametrized FP action 

The following table collects the numerical values of the non-linear and linear 
parameters describing the approximate FP action in the range of lattice spacing 
0.03 fm < a < 0.3 fm. 

The set of parameters consists of four non-linear parameters 77^*^^ , cf^ , c^'^ , 03"' 
describing the asymmetrically smeared links wjf'^ and fourteen linear param- 
eters pki with < fc + ^ < 4. This set approximates reasonably well the true 
FP action in the range of a given above. (For smaller fluctuations - in the 
intermediate steps of the parametrization - we used polynomials for ri{x) and 
Ci{x) up to forth order.) 

The optimal non-linear parameters are found to be 

77(°) = -0.038445 , cf ^ = 0.290643 , 4"' = -0.201505 , 4"' = 0.084679 . 



The linear parameters are collected in table A.l 
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/ = 


I = 1 


/ = 2 


I = 3 


/ = 4 


k ^ 







0.442827 


0.628828 


-0.677790 


0.176159 


k = 


1 


0.051944 


-0.918625 


1.064711 


-0.275300 




k = 


2 


0.864881 


-0.614357 


0.165320 






fc = 


3 


-0.094366 


-0.020693 








k = 


4 


0.022283 











Table A.l: Linear parameters pki of the parametrized FP action. 



B Variational techniques 

In a Monte Carlo simulation we measure the N x N correlation matrix 

Co.i3{t) ^ {0\Oo.{t)Ol{Q)\0) . (B.l) 

To determine the coefficients Va of the linear combination ^aOa which 

has the largest overlap to the ground state relative to the excited states one has 
to minimize the effective mass given by 



m{to,ti 



In 



{v,C{h)v) 



/{ti~to). 



(B.2) 



[iv,C{ta)v) 

The vector v is obtained by solving the generalized eigenvalue equation p6| , |67| 

C{h)v = Xito,h)C{to)v , (B.3) 

where < to < ^i- 

Assume ffist that only the lowest lying N states contribute to C{t), i.e. 



Ca/3 (t) 



N 



(B.4) 



where Ei < E2 < ... < En arc the energy levels in the given symmetry 
channel and tpna — {0\Oa\n) is the "wave function" of the corresponding state. 
The solution of eq. (B.3) is given by the set of vectors {vn} dual to the wave 
functions, i.e. {vn^ipm) = Snm- Multiplying eq. (B.4) by Vn one obtains 



C{t)vn = e 



E,,t, 



= e-^"(*-*«)C(fo)t'n. (B.5) 



This gives A„(io,ii) — exp(— £'„(ti — to)) for the eigenvalues in eq. (B.3). Of 
course, contributions from states with n > N and statistical fluctuations distort 



eq. (B.4), therefore the stability of eq. (B.3) is an important issue. 

Observe that eq. (B.3) is well defined only for positive definite C{to). Because 
of statistical fluctuations, however, the measured correlation matrix C(io) is not 
necessarily positive for to > 0. This is the reason why one usually considers 
only the to = case in applying the variational method, especially with a large 
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number of operators. On the other hand, it is obvious that C(0) is contaminated 
by highly excited states and contains only restricted information on the low 
lying part of the spectrum. Therefore it is desirable to take to > 0. This can be 
achieved in the following way [ |68[ . 
We first diagonalize C(to), 

C{to)ip; = \^^^, Xi>...>Xn, (B.6) 

and project the correlation matrices to the space of eigenvectors corresponding 
to the M highest eigenvalues, 

C^'{t)^{^,,C{t)^,), t,j = l,...,M. (B.7) 

By choosing the operator space too large we introduce numerical instabilities 
caused by very small (even negative) eigenvalues with large statistical errors 
due to the fact that the chosen operator basis is not sufhciently independent 
on the given MC sample. By choosing M appropriately we can get rid of those 
unstable modes while still keeping all the physical information. In this way we 
render the generalized eigenvalue problem well defined. 

Of course the final result should not depend on the choice of M and one 
has to take care in each case that this is really the case. Our observation is 
that for any acceptable statistics one always finds a plateau in M for which the 
extracted masses are stable under variation of M. 

In a next step we determine the vectors Vn, n = 1,... ,M through the 
generalized eigenvalue equation in the truncated basis: 

C^HhW = e-^"(*i-*°)C^^(to)z^„ . (B.8) 
This equation yields the spectrum i?„. However, the procedure - although it is 



exact for a correlation matrix which has exactly the form in eq. (B.4) ~ is highly 
non-linear, and a small statistical fluctuation can be enhanced by it and cause 
a systematic shift in the energy values obtained, even when the instabilities are 
avoided by the truncation to M < N. 

In order to avoid this pitfall we use the (approximate) dual vectors u„ ob- 



tained from eq. (B.J) to restrict the problem to a smaller, therefore more stable 
subspace. 

Define the new correlation matrix of size K x K (with K < M) hy 

C^{t)^iv,,C^'{t)v,), ,K<M. (B.9) 

The steps performed until now can be thought of as a preparation for choosing 
the appropriate set of operators, i.e. linear combinations of original Oa operators 
which effectively eliminate the higher states. The correlation matrix Cfj (t) is 
then considered as a primary, unbiased object. 

The next step is to fit Cfj (t) in the range t = tmin ■ ■ ■ ^max using the ansatz 

K 

Cf^ it- {V, E}) = e-'^-VmC, , (B.IO) 
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where tpni, En are the free parameters to be fitted. 

Usually we choose K = 1 and 2. For the Af^ glueball, however, K = 2 and 
3 are chosen since we do not subtract the vacuum contribution {Oa){Op)* from 
the correlators but consider instead the vacuum state together with the glueball 
states in this channel (cf. remarks in section S.3). 

In the fitting step we use a correlated fit which takes into account the cor- 
relation between (t) and C^j,{t'), i.e. using the inverse of the corresponding 
covariance matrix Gov{i, j,t;i' , j' ,t') as a weight in the definition of x^. This 
has the advantage over the usual (uncorrelated) x^ that the value of the latter 
can be artificially small if the quantities to be fitted are strongly correlated. 
Note however, that (as usually with sophisticated methods) the correlated x^ 
fit can have its own instabilities if the number of data is not sufficiently large 



C Simulation parameters and results 
C.l Deconfining phase transition 



lattice size 


P 


sweeps 


Tp 


Tint 


2 X 10^ 


2.3550 


30000 




260 




2.3575 


30000 


4300 


283 




2.3560 


30000 


4600 


280 


2x8^ 


2.3300 


14240 




29 




2.3500 


10144 




93 




2.3550 


5120 




127 




2.3575 


12288 


1400 


202 




2.3700 


10144 




114 


2x6^ 


2.3250 


8096 




35 




2.3500 


14144 


650 


105 




2.3600 


10000 


700 


96 




2.3750 


10144 
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Table C.l: Run parameters of the critical temperature simulations at Nr = 2 
including the persistence time Tp and the integrated autocorrelation time Tint of 
the Polyakov loop operator. 
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lattice size 




sweeps 


Tp 


Tint 


3 X 12^ 


2.675 


25000 




114 




2.680 


45000 


3200 


188 




2.685 


24000 




96 




2.690 


20000 




53 


3 X 10^ 


2.670 


18000 




67 




2.680 


42000 


2300 


89 




2.685 


48000 


2400 


104 




2.690 


27000 




85 


3x8^ 


2.650 


10096 




43 




2.660 


10000 




48 




2.670 


26000 




41 




2.680 


30000 


1400 


64 




2.690 


19000 




53 




2.710 


10000 




35 



Table C.2: Run parameters of the critical temperature simulations at A'^^ = 3 
including the persistence time Tp and the integrated autocorrelation time Tint of 
the Polyakov loop operator. 



lattice size 


P 


sweeps 




Tint 


4 X 14^ 


2.917 


50405 


4300 


62 




2.922 


51812 


4700 


67 




2.930 


44607 




64 


4 X 12=* 


2.850 


15000 




19 




2.890 


15000 




30 




2.910 


33000 




34 




2.920 


33000 


3700 


66 




2.930 


15000 




38 


4 X 10=* 


2.850 


10000 




22 




2.880 


16000 




37 




2.890 


21124 




18 




2.900 


35000 




34 




2.910 


35000 


2100 


36 




2.920 


20000 
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Tabic C.3: Rim parameters of the critical temperature simulations at TV,- = 4 
including the persistence time Tp and the integrated autocorrelation time Tint of 
the Polyakov loop operator. 
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C.2 Static quark-antiquark potential 



(3 lattice volume lattice size [fm] # measurements 

TIoo T¥ E45 43 X 90 

3.150 12^ 1.61 42 X 50 

2.927 14^ 2.39 40 x 40 

2.860 10^ 1.84 43 x 90 

2.680 12^ 2.72 51 x 40 

2.361 12^ 4.02 57 x 40 



Table C.4: Run parameters for the simulations of the static quark-antiquark 
potential. 





fit range 


aVo 


a 






3.400 


2 - 


6 


0.781(1) 


0.251(9) 


0.063(1) 


1.02 


3.150 


2 - 


5 


0.820(15) 


0.285(20) 


0.099(3) 


0.75 


2.927 


2 - 


6 


0.812(16) 


0.272(20) 


0.161(3) 


1.35 


2.860 


1 - 


4 


0.801(5) 


0.262(3) 


0.189(2) 


1.17 


2.680 


1 - 


4 


0.777(5) 


0.255(4) 


0.287(2) 


0.43 


2.680 


2 - 


6 


0.778(41) 


0.256(54) 


0.287(7) 


0.65 


2.361 


1 - 


4 


0.615(11) 


0.179(8) 


0.629(4) 


0.99 


2.361 


2 - 


5 


0.59(11) 


0.15(13) 


0.634(22) 


1.41 



Table C.5: Results from global correlated fits of the form (|o|) to the static 
quark potentials. The second column indicates the fit range in r and the last 
column per degree of freedom, x^/A^df- 
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P 


r 




fit range 


V{r) 




O A A A 

o.4(J(J 


1 


5 


1 - 





0.5874(2) 


0.76 




o 

/ 


r 






D 


n TO A /I c \ 

U.7oU4(oJ 


Z.i9 




Q 

o 


c: 



Q 

o - 





A QQc:/o^ 
U.ooO^Z ) 


1 07 
l.Z / 




/I 
4 


Q 

o 


q 

o - 





A A/?A/Q\ 

U.yoy^o J 


i.Z4 




r 
t) 


4 





D 


i.U4D(4J 


u.yi 




D 


4 


Z - 


rr 



i.liO^O j 


U.oo 




/ 


Q 

o 


6 - 





1.17(2j 


U.io 


o.ioU 


i 


r 




3 - 


5 


A C /I A C f 0\ 

U.d4Uo(oJ 


0.77 




2 


4 


z - 


D 


0.8756(5) 


0.63 




Q 
O 


4 


Z - 





i.Uzo^i j 


U.D4 




4 


Q 
O 








i.i4^(zj 


U.io 




tr 



Q 

o 








i.ZOO^O j 


U.o4 




D 


Q 

o 




z - 





1 QC ^ 1 ^ 


i.Uo 


z.yzY 


i 


4 


o 


( 


U.7Uoz(zj 


0.42 




o 
z 


Q 
O 


o 

z - 


7 




U.DO 




3 


3 


z - 


7 


1 AAA / A\ 

1.202(2) 


0.56 




4 


4 




z - 





1.383(5) 


A O 1 

0.31 




c: 



Q 
O 


o 

z - 


7 
1 


i.0DU(O j 


U.oi 




b 


Q 
O 


Z - 


r 




1.71(2) 


O.oZ 




7 


o 


9 

z - 


a 

D 


i.yz(3j 


l.Zo 


Z.oDU 


1 


Q 
O 


o 

z - 


A 

4 


U. / ZD / (4j 


1 c;n 






o 
O 


i - 


4 


1.047(1) 


0.56 




Q 
C) 


4 


i - 


4 


i.z ^o(z ) 


U.oo 




A 

4 


o 
I 


Z - 


/I 
4 


1.4(5(5(0 ) 


U.oU 




r- 




6 


Z - 


A 

4 


1.67(z) 


U.6o 


Z.DOU 


1 

i 


A 

4 


9 

z - 


D 


0.8091(0) 


U.zi 







A 

4 


2 - 


6 


i.zzony ) 


u.yo 




3 


4 


2 - 


6 


1.553(3) 


0.33 




4 


3 


1 - 


5 


1.862(3) 


0.33 




5 


2 


2 - 


6 


2.15(3) 


0.89 




6 


2 


2 - 


5 


2.51(8) 


0.14 


2.361 


1 


3 


2 - 


5 


1.0641(6) 


0.33 




2 


3 


1 - 


6 


1.783(1) 


0.31 




3 


2 


1 - 


5 


2.443(4) 


0.75 




4 


2 


1 - 


6 


3.09(2) 


0.84 




5 


1 


1 - 


5 


3.73(6) 


2.31 




6 


1 


1 - 


6 


4.5(3) 


0.44 



Table C.6: Potential values extracted from fits of the form Z[r) exp(— ty(r)) to 
the ground state of the string correlators. For each /?-value and distance r we 
list the plateau regions (fit range tmin — tmax), the extracted potential values 
y(r) and the per degree of freedom, x^/-^df- Note that = \ and t\=2 
was chosen in all cases. The column entitled with M denotes the number of 
operators kept after the first truncation. 
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C.3 Glueballs 



/3 


lattice 


a[fm] 


ro/a 


# measurements 


3.40 


14^ 


0.10 


4.83(4) 


206 X 70 


3.15 


12^ 


0.13 


3.72(2) 


202 X 50 


2.86 


10^ 


0.18 


2.74(1) 


160 X 50 



Table C.7: Run parameters of the glueball simulations. Values for the cou- 
pling /3, the lattice size and the obtained statistics are listed. The estimate of 
the hadronic scale ro in terms of the lattice spacing a is given as well as the 
approximate lattice spacing in fermi. 



Channel 


to/ti 


M 


fit range 


X7^DF 


energies 




1/2 


3 


1 - 4 


0.02 


1.41(10) 








2 - 4 


0.02 


1.40(38) 




0/1 


25 


1 - 4 


0.56 


1.38(8) 








2 - 4 


0.36 


1.50(40) 



Table C.8: Results from fits to the (3 = 2.86 glueball correlator on the 10^ 
lattice. Only five loop shapes were measured on 5 different smearing schemes. 
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Channel 


to/h 


M 


fit range 


XV^DF 


energies 


Ar 


1/2 


5 


1 - 3 


0.61 


1.03(3) 








2 - 3 


0.00 


1.10(10) 








1 - 4 


2.02 


1.03(3) 








2 - 4 


1.07 


1.12(11) 




0/1 


25 


1 - 4 


1.62 


1.02(3) 








2 - 4 


0.02 


1.12(9) 


E++ 


1/2 


4 


1 - 3 


1.26 


1.53(6) 




0/1 


48 


1 - 3 


1.41 


1.46(5) 


T++ 


1/2 


4 


1 - 3 


0.68 


1.61(6) 




0/1 


48 


2 - 4 


1.32 


1.83(23) 




1/2 


3 


1 - 3 


0.84 


1.65(18) 


E-+ 


1/2 


3 


1 - 3 


0.00 


1.97(20) 




0/1 


15 


1 - 3 


0.09 


2.06(16) 




1/2 


5 


1 - 3 


0.00 


1.39(27) 




0/1 


22 


1 - 3 


0.00 


1.92(11) 




1/2 


4 


1 - 3 


2.70 


2.10(18) 




0/1 


25 


1 - 3 


0.05 


2.04(12) 



Table C.9: Results from fits to the = 3.15 glueball correlators on the 12^ 

lattice: to/ti arc used in the generalized eigenvalue problem, M denotes the 
number of operators kept after the truncation in C{to)- 
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Channel to/ti M fit range X/^^df energies 



A++ 
^1 


1 /o 


fi. 




1 
i 


A 

- 4 


u. / y 






n/1 

U/i 




1 

± 


A 






E++ 


1/2 


11 


1 


- 4 


0.03 


1.23(5) 






8 


1 


- 4 


0.19 


1.27(3) 




0/1 


60 


1 


- 4 


0.02 


1.23(2) 


T++ 


1/2 


5 


1 


- 4 


0.40 


1.23(3) 






7 


1 


- 4 


0.16 


1.20(3) 




0/1 


48 


1 


- 4 


1.16 


1.25(2) 




1/2 


3 


1 


- 3 


0.24 


1.40(9) 




0/1 


15 


1 


- 3 


0.12 


1.46(5) 






15 


2 


- 4 


0.10 


1.38(20) 


E-+ 


1/2 


3 


1 


- 3 


0.34 


1.68(7) 




1/2 


4 


1 


- 3 


0.09 


1.63(7) 




1/2 


8 


1 


- 3 


2.49 


1.64(16) 






6 


1 


- 3 


0.17 


1.76(10) 




0/1 


25 


1 


- 3 


0.07 


1.65(6) 



Table C.IO: Results from fits to the = 3.40 glueball correlators on the 14^ 
lattice: t^jtx are used in the generalized eigenvalue problem, M denotes the 
number of operators kept after the truncation in C(<o)- 





/3 = 2.86 


/3 = 3.15 


/3 = 3.40 




1.41(10) 


1.05(6) 


0.84(2) 






1.53(6) 


1.23(5) 






1.61(6) 


1.23(3) 






1.65(18) 


1.40(9) 






1.97(20) 


1.68(7) 






1.92(11) 


1.63(7) 


T+- 




2.10(18) 


1.64(16) 



Table C.ll: Final glueball mass estimates in terms of the lattice spacing, ama- 
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J ,J = 2.8() 


J = 3.1o 


J = 3.40 




3.87(27) 


3.92(23) 


4.04(12) 


E++ 


2 


5.70(23) 


5.96(24) 


T++ 


2 


5.98(21) 


5.96(14) 







6.13(67) 


6.74(42) 


E-+ 


2 


7.32(74) 


8.12(35) 




2 


7.14(41) 


7.88(35) 


T+- 


1 


7.81(67) 


7.93(78) 



Table C.12: Final glueball mass estimates in terms of ro, romo- The continuum 
spin interpretation of each channel is denoted by J. 
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